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Abstract 

We reanalyze critically the generalized factorization hypothesis in non-leptonic two- 
body B-decays discussed recently by several authors. In particular we address the 
determination of the factorization scale \Xf and of the non-perturbative parameters 
^ NF (mb) and £,2 F ( m b) which are supposed to measure non-factorizable contributions 
to hadronic matrix elements with £ 4 NF (/i/) = 0. We emphasize that both u,f and 
^ NF (mb) are renormalization scheme dependent and we demonstrate analytically and 
numerically that for any chosen scale \if — 0(m^) it is possible to find a renormal- 
ization scheme for which ^ F (/i/) = Cf^A*/) — 0- The existing data indicate that 
such "factorization schemes" differ from the commonly used schemes NDR and HV. 
Similarly we point out that the recent extractions of the effective number of colours 
N from two-body non-leptonic B-decays while /i and renormalization scheme inde- 
pendent suffer from gauge dependences and infrared dependences. 
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1 Introduction 



Two-body non-leptonic B-decays play an important role in the phenomenology of 
weak decays not only probing the structure of weak interactions corrected by short 
distance QCD effects but also providing some insight into the non-perturbative phe- 
nomena related to long distances. The increasing experimental information on these 
decays, in particular from the CLEO detector, stimulated recently several new theo- 
retical analyses of these decays. The most extensive analyses of this type are based 
on the factorization of hadronic matrix elements of local operators |]]]-[f§] which has 
recently been extended to the so-called generalized factorization hypothesis @j-[P]. 

In the strict factorization approach two-body decays are parametrized in terms of 
two phenomenological parameters a\ and a 2 which in QCD are given by 

ai(//) = Ci(//) + jjC 2 (n) , a 2 (» = C 2 (ji) + jjC x (ji). (1) 

Here C\,i{yi) are the short distance Wilson coefficient functions of the relevant current- 
current operators Ox^i for which explicit expressions will be given below. N is the 
number of colours with N = 3 in QCD. 

One distinguishes then three classes of decays for which the amplitudes have the 
following general structure f|. |J : 

^4i = ^WiM(Oi}f (Class I) (2) 

Au = ^V CKM a 2 (fJL)(0 2 ) F (Class II) (3) 
G F 

Am. = —^VckmWM +xa 2 (iJ.)](0 1 ) F (Class III) (4) 

Here Vckm denotes symbolically the CKM factor characteristic for a given decay. 
(Oi) F are factorized hadronic matrix elements of the operators 0« given as products 
of matrix elements of quark currents and x is a non-perturbative factor equal to unity 
in the flavour symmetry limit. 

The simplicity of this approach is very appealing. Once the matrix elements 
(Oi) F have been expressed in terms of various meson decay constants and generally 
model dependent form factors, predictions for non-leptonic heavy meson decays can 
be made. An incomplete list of analyses of this type is given in [|, |IU| and will be 
extended below. 
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On the other hand, it is well known that non-factorizable contributions must be 
present in the hadronic matrix elements of the current-current operators 0± and O2 
in order to cancel the /1 dependence of Cj(/i) or aj(/i) so that the physical amplitudes 
do not depend on the arbitrary renormalization scale \x. (Oi)p being products of 
matrix elements of conserved currents are /^-independent and the cancellation of the 
\i dependence in does not take place. Consequently from the point of view 

of QCD the factorization approach can be at best correct at a single value of /1, the 
so-called factorization scale fif. Although the approach itself does not provide the 
value of /i/, the proponents of factorization expect /i/ = 0(mb) and /i/ = 0(m c ) for 
B-decays and D-decays respectively. 

The fact that (Oi)p are /i-independent but are /i-dependent, which is clearly 
inconsistent, inspired a number of authors [f|-|| to generalize the concept of factor- 
ization. 

In the formulation due to Neubert and Stech |7| the /i-dependent parameters 
ai(/i) and a 2 (/t) are replaced by /i-independent effective parameters af ff and a| ff . The 
latter depend on Cj(/t) and two non-perturbative parameters £i(/i) and Eg([i) which 
parametrize the non-factorizable contributions to the hadronic matrix elements of the 
operators 0\p- In the case of strict factorization e$ vanish and a\ s 2 reduce to ai^/i). 
The /t dependence of £i(/i) cancels the /i-dependence of Ci )2 (/i) so that af^ are indeed 
scale independent. 

From the phenomenological point of view there is no change here relative to the 
standard factorization as only a,(/i) have been replaced by af in the formulae (^])-(f|). 
On the other hand, as stressed in J7[, the new formulation should allow in principle 
some insight into the importance of non-factorizable contributions to hadronic matrix 
elements. 

In this context we should remark that in the recent literature mainly the \x- 
dependence of the non-factorizable contributions has been emphasized. Their renor- 
malization scheme dependence has often not been discussed. It is the latter issue 
which will be important in the discussion below. Indeed at the next-to-leading level 
in the renormalization group improved perturbation theory the coefficients Cj(/i) 
depend on the renormalization scheme for operators. Again only the presence of 
non-factorizable scheme dependent contributions in (0») can remove this scheme de- 
pendence in the physical amplitudes and in particular in a\ . The renormalization 
scheme dependence emphasized here, and discussed in the context of strict factoriza- 
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tion in JTT], is rather annoying from the factorization point of view as it precludes 
a unique phenomenological determination of /i/ as we will show explicitly below. In 
particular we will demonstrate that for any chosen scale fif = 0(m^,) it is possible 
to find a renormalization scheme for which the non-factorizable parameters £i i g(/xj) 
simultaneously vanish. This finding casts some doubts on the usefulness of the formu- 
lation in |7j with respect to the study of non-factorizable contributions to non-leptonic 
decays. 

The generalized factorization presented in || ^ || is similar in spirit but includes 
more dynamics than the formulation in M . Here the non-factorizable contributions to 
the matrix elements are calculated in a perturbative framework at the one-loop level. 
Subsequently these non-factorizable contributions are combined with the coefficients 
Ci(n) to obtain effective fi and renormalization scheme independent coefficients C^ s . 
The effective parameters a? are given in this formulation as follows: 

with analogous expressions for a| ff (i = 3 — 10) parametrizing penguin contribu- 
tions. Here N eS is treated as a phenomenological parameter which models the non- 
factorizable contributions to the hadronic matrix elements. In particular it has been 
suggested in |5|, [8], [J that the values for N eS extracted from the data on two-body 
non-leptonic decays should teach us about the pattern of non-factorizable contribu- 
tions. 

Unfortunately, as we will demonstrate below, also this approach has its weak 
points. Although Cf s 2 are A* an d renormalization scheme independent, they are both 
gauge and infrared regulator dependent. The latter dependences originate in the 
perturbative evaluation of the scheme dependent finite contributions to the matrix 
elements, needed for the cancellation of the renormalization scheme dependence of 
Ci(fi). Consequently, whereas the extracted N is renormalization scheme and renor- 
malization scale independent, it is a gauge and infrared regulator dependent quantity. 
This finding casts some doubts on the usefulness of the formulation in || H |9|] with 
respect to the study of non-factorizable contributions to non-leptonic decays. 

The rest of our paper amounts to putting all these statements in explicit terms. 
In section 2 we review the approach in 0. In section 3 we reformulate this approach 
by replacing the parameters Ei^in) by two new parameters £ff . This reformulation 
allows us to make our points with regard to in a more transparent manner. In 
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particular we derive general expressions which allow to find, for a given [if, the 
renormalization scheme in which £i^{fi) or £f f (/i) simultaneously vanish. In section 
4 we illustrate our points with a few numerical examples and in section 5 we make a 
critical analysis of the approach in [|5], [8], [|. We end our paper with a brief summary 
and conclusions. 



2 Generalized Factorization 

In order to describe generalized factorization in explicit terms let us consider the 
decay B° — > D + n~ . Then the relevant effective Hamiltonian is given by 

H eS = ^W&ICMaOOx + cmo 2 ] , (6) 

where 

Oi = (d a u a ) v _ A (cpbp) v _ A , 2 = (d a up)v-A(cf3b a ) v _ A (7) 

with (a, (3 = 1, 2, 3) denoting colour indices and V — A referring to 7^(1 — 75). Ci(/i) 
and C2(/i) are short distance Wilson coefficients computed at the renormalization 
scale ji = 0{m b ). Note that we use here the labelling of the operators as given in 0. f|] 
which differs from ||12|- ||15|| by the interchange 1^2. Since all four quark flavours 
entering the operators in (|7|) are different from each other, no penguin operators 
contribute to this decay. 

Using Fierz reordering and colour identities one can rewrite the amplitude for 
B° D+tt- as 

A(B° -> D+O = ^W&af (O^p , (8) 

where 

(Oi) F = (tt- I (du) V -A I 0>(£> + I (c6)y_A I B°) (9) 

is the factorized matrix element of the operator 0\ and summation over colour indices 
in each current is understood. 

The effective parameter a\ s is given by |7J 

af = (Ci(/i) + ^(a«)) [1 + e^\n)\ + C 2 (M BD ^(V)- (10) 
£1 D (/i) and Eg BD ' n \ij) are two hadronic parameters defined by M 

(Oi)f 



5 



and 

(BD,n), v _ D + \(dt a u) V -A(ctab) V - A \B ) 

£ 8 = 1 TTTT \ l l 

with i a denoting the colour matrices in the standard Feynman rules. £j(/x) parametrize 
the non-factorizable contributions to the hadronic matrix elements of operators. In 
the case of strict factorization £j vanish. 

It should be emphasized that no approximation has been made in (Q). Since the 
matrix element (0\)f is scale and renormalization scheme independent this must also 
be the case for the effective coefficient af ff . Indeed the scale and scheme dependences 
of the coefficients Ci(/x) and C^/x) are cancelled by those present in the hadronic 
parameters £j(/x). We will give explicit formulae for these dependences below. 

A similar exercise with the amplitude for B° — > D°tt° gives 

A(B° -> D%°) = ^=V cb V: d af(0 2 ) F , (13) 

where 

(0 2 ) F = (D° | (cu) v „ A | 0)(7T° I {db) V -A I B°) (14) 

is the factorized matrix element of the operator Oi- 
The effective parameter a| ff is given by pj 



af 



(CM + ^CxC/i)) [1 + ei^ifi)) + C x {^t' D \n) . (15) 
£ ^ Bw ' D ) ^ anc [ £ ^> D ) are two hadronic parameters defined by 

,r D V) - &^mi^m _ , (16) 

{O2/F 

and 

(fl ^) (7r° J D |(a a «) y _ A (Ja)y-A|S°) 

4 W = 2 77T~\ (17) 

Again the /x and scheme dependences of £j in (|16|) and (|17 ) cancel the corresponding 



dependences in Cj(/x) so that the effective coefficient a| ff is /x and scheme independent. 



Following section 5.1 of |j2j it is straightforward to find the explicit /x and scheme 
dependences of the hadronic parameters £,(/x). To this end we note that the /x de- 
pendence of the matrix elements of the operators 0± = {0% ± 0-i)j1 is given by 

H 

(0±(/x)) = tf±(m b ,/x)(0 ± (m b )) (18) 
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where the evolution function £/±(m b , fi) including NLO QCD corrections is given by 



U±(m h ,fi) 



a s (m h ) 

1 + ^r J± 







d± - 




_a s (m h )_ 





An 



-J4 



with 



d 



± 



7 i 0) = ±2(3 T l) 



7i 1} 



Pa 

3=Fl 
6 



11-?/ 

3 J 



(0) 
7± 

2/3o ' 

Pi = 102 - f / 



-21 ± -/ - 2ftK± 



(19) 

(20) 

(21) 
(22) 



Here k±, introduced in ||11|| , distinguishes between various renormalization schemes: 

(23) 



(NDR) [16] 
(HV) 



Thus J± in 



^6-3 (DRED) JT7 
can also be written as 

3=Fl 



(o) 
7± 



J± — (J±)ndr H 7 — K ± — (J±)ndr i -ttt^± 



6 
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The MS 1 coupling [IS| is given by 

An 

a s (fjL) 



1 - 



Pi lnln(/i 2 /A|^ ; 



/3 2 ln(/i 2 /A 



A/5 7 



(24) 



(25) 



The formulae given above depend on /, the number of active flavours. In the case of 

(26) 



B-decays / = 5. The present world average for a s (Mz) is [ |I~9|j : 

a s (M z ) = 0.118 ± 0.003 = (225 ± 40) MeV 



where the superscript stands for / = 5. 

Having these formulae at hand it is straightforward to show that the //-dependence 
of Ei(ft) and £s(/f) is governed by the following equations: 

1 



l + £l(A0 



+ 



1 + 
1 



N 
1 

N 



[1 +£i(m b )] +e 8 (m h ) 
[1 +ei(m b )] -e 8 (m b ) 



U + (m hl n) 
f/_(m b ,/i) 



(27) 



1 

2 
1 

+ 2 



1 • -8(m b )-(l- — ) I 
1 • y) £8(m b )-(l- — ) I 



U + (m h ,n) 
U-(m h ,(j) 



(2f 



These formulae reduce to the ones given in |7J when J± in flT9"| ) are set to zero. 
They give both the //-dependence and renormalization scheme dependence of £j. The 
latter dependence has not been considered in [^]. 



3 A Different Formulation 



In order to be able to discuss the relation of our work to the one of JF| we have used un- 
til now, as in 0], the hadronic parameters Si(fi) and Es(fj) to describe non-factorizable 
contributions. On the other hand, it appears to us that it is more convenient to work 
instead with two other parameters defined simply by 



af 



a 2 (/i)+£ 2 NF 0u) 



(29) 



where ai(/i) and a 2 (fi) are given in ([!]). 
Comparison with fllCf) and (O) gives 



(30) 
(31) 



where 



e 



{Btv,D) 



(32) 
(33) 



£i(//) = e\ ' e 8 (//) 

eiW = e i > £ s(w 

and ai(n), given in ([]]), are the parameters used in the framework of the strict fac- 
torization hypothesis in which £^ F (/x) are set to zero. Their /i and renormalization 
scheme dependence can be studied using 

z+(ji) + z-(p) 



where 



with 



Ci(m) 



z±(n) 



1 + 



2 

47T 



cm 



z + {n) - z-(fi) 





~a s (M w )~ 


d± - 









Air 



fl± = ^[±ll + /c±] 
o 



(34) 

(35) 
(36) 



and all other quantities defined before. The /i and scheme dependences of £f can in 
principle be found by using the dependences of C;(/^) given above and Ei(fi) in $27\) 
and (|28|). To this end, however, one needs the determination of the non-perturbative 
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parameters and £«(//) at a single value of /i. If, as done in |7]], af are assumed 
to be universal parameters, the determination of £i(/i) and £j(/i) is only possible if 
one also makes the following universality assumptions: 

ei(ii) = e 8 (/x) = e 8 (/i) . (37) 

In J7| such an assumption was not necessary since £i(/x), £i(/f) and £s(aO were ne- 
glected, and only e 8 (fi) was kept in the analysis. 

With the assumptions in (|37|), £i(/x) and £s(/-*) can indeed be found once the 
effective parameters af have been determined experimentally. Using (|l(J) and (|i~5l) 
together with ( |37l) we find 

CMaf - C 2 (ji)a? 

£l(/i) = ^F^) — (38) 

^-^-(SS^)^^ (39) 

On the other hand £f F (/x) can be determined without the universality assumption 
( p7[ ) from two decays simply as follows 

ef F (^) = of - aiM e 2 NF (/i) = af - a 2 (//) (40) 



We will analyze the formulae (|3"8D, ( |59]) and ( [TOP in the next section. 

The formulae in ( |40"1) make it clear that the strict factorization in which ^ NF (yu) 
vanish can be at best correct at a single value of fi, the so-called factorization scale 
fif. In the first studies of factorization fif = has been assumed. It has been 
concluded that such a choice is not in accord with the data || [TD|. |j] . 



The idea of the generalized factorization as formulated by Neubert and Stech J7J 
(see also || || for earlier presentations) is to allow Hf to be different from m b and to 
extract first the non-factorizable parameters £j(m b ) from the data. Subsequently the 
factorization scale /// can be found by requiring these parameters to vanish. 

In the numerical analysis of this procedure done in [0 a further assumption has 
been made. Using large N arguments it has been argued that £i(p) can be set to 
zero while £s(aO can be sizable. The resulting expressions for af ff are then 



- Ci(m b ), af = a 2 (m h ) + Ci(m b )e 8 (m b ) (41) 

where additional small terms have been dropped in order to obtain the formula for 
af. Using subsequently the extracted value af = 0.21 ± 0.05 together with the 
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coefficients Ci(m\>) from |TjJ one finds £s( m b) = 0.12 ± 0.05 Next assuming 
£8 (a*/) = one can find the factorization scale /i/ by inverting the formula 

4a s (m b ) m b 
£s ("lb) = ^ In (42) 

which follows from (|28| ) with s 8 (fif) = and £i(m b ) = 0. Thus 

37T£ 8 (m b ) 



/i/ = m b exp 



4a s (m b ) 



(43) 



Taking m b = 4.8 GeV and a s {m^) = 0.21 (corresponding to a s {Mi) = 0.118) we find 
using e 8 ( m b) = 0.12 ± 0.05 a rather large factorization scale /i/ = (15.9^^) GeV, 
roughly a factor of 3-4 higher than m b . This implies that non-factorizable contri- 
butions in hadronic matrix elements at scales close to m b are sizable. This is also 
signaled by the value of Eg(mb) ~ 0.12 which is as large as the factorizable contribu- 
tion a2{m\ ) ) = 0.10 to the effective parameter a| ff = 0.21 ± 0.05. 

We would like to emphasize that such an interpretation of the analysis of Neubert 
and Stech would be misleading. As stressed in JTT[] the coefficient a<i (/i) is very 



strongly dependent on the renormalization scheme. Consequently for a given value 
of a| also £,2 F ( m h) an d e^{m^) are strongly scheme dependent. This shows |~1 



that a meaningful analysis of the /^-dependences in non-leptonic decays, such as the 
search for the factorization scale fif, cannot be be made without simultaneously 
considering the scheme dependence. This is evident if one recalls that any variation 
of \ij in the leading logarithm is equivalent to a shift in constant non-logarithmic 
terms. The latter represent NLO contributions in the renormalization group improved 
perturbation theory and must be included for a meaningful extraction of fif or any 
other scale like A^r. However, once the NLO contributions are taken into account, 
the renormalization scheme dependence enters the analysis and consequently the 
factorization scale fif at which the non-factorizable hadronic parameters or 
Ei(nf) vanish is renormalization scheme dependent. 

From this discussion it becomes clear that for any chosen scale /x/ = (9(m b ), it is 
always possible to find a renormalization scheme for which 

t? F <J*f) = ^2 F M = • (44) 

Indeed as seen in (f40"D ^ F (/i) depend through cij(/z) on k± which characterize a 
given renormalization scheme. The choice of k± corresponds to a particular finite 
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renormalization of the operators 0± in addition to the renormalization in the NDR 
scheme. It is then straightforward to find the values of k± which assure that for a 
chosen scale /i/ the conditions in ([14]) are satisfied. We find 



K+ 



cff 



3 af + af 

4 w + m 



- 1 



An 



3 a? 



cff 



,cff 



2 [2 W-(ji f ) 



An 



3(J 4 



a s {ii f ) 2 



J NDR 



WDR 



where 



~a s (M w )~ 


d± - 







l+ MMWl {B± _ J±) 

Att 



(45) 
(46) 

(47) 



with (J±)ndr being the values of J± in the NDR scheme. W±(jif) are clearly renor- 
malization scheme independent as B± — J± are scheme independent. 



4 Numerical Analysis 



Before presenting the numerical analysis of the formulae derived in the preceding 
section, we would like to clarify the difference between the Wilson coefficients in ( [34]) 
and (|3~5|) used by us and the ones employed in 0. In |7| the scheme independent 
coefficients z±([i) of |16[] instead of z±(fj) have been used. These are obtained by 
multiplying z±(fj) by (1 — B±a s (/i) / An) so that 



a s (M ] 



w , 



1 + 



a s (M w ) - ajjj) 
An 



(B± - J±) 



(48) 



These coefficients are clearly not the coefficients of the operators 0± used in and 
here. In order to be consistent the matrix elements (0±) should then be replaced by 



(d ± ) = (l + B ± a s (ii)/An)(0 ± ). 



(49) 



This explains why the results of our numerical analysis differ considerably from the 
ones presented in @. We strongly advice the practitioners of non-leptonic decays not 



to use the scheme independent coefficients of |T6]] in phenomenological applications. 
These coefficients have been introduced to test the compatibility of different renor- 
malization schemes and can only be used for phenomenology together with (0±). 
This would however unnecessarily complicate the analysis and it is therefore advis- 
able to work with the true coefficients Cj(/i) of the operators Oi as given in (]34[) and 
351). 
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Figure 1: ei t s(fi) in the NDR and HV schemes. 

In the following values of af 1 have been extracted from existing data on two- 
body B-decays 

af = 1.08 ± 0.04 af = 0.21 ± 0.05 (50) 

with similar results given in [5|, || [TjJ ||, [| . In order to illustrate various points made 
in the preceding section, we take the central values of a| ff in ([50]). Using (|38| ) -(|4*0| ) we 
calculate £j(/i) and as a function of /i in the range 2.5 GeV < fi < 10 GeV 

for the NDR and HV schemes. The results are shown in fig. [I] and fig. |^. We observe 
that and £f F (/i) are only weakly \i and scheme dependent in accordance with 

the findings in [jllj, where these dependences have been studied for aj(/x) defined in 
([I]). The strong [i and scheme dependences of a^G") found there translate into similar 
strong dependences of e%{\i) and £^ F (/i). 
We make the following observations: 

• £i(/u) and £j^ F (/i) are non-zero in the full range of /i considered. 

• es(fi) and £^ F (/i) vary strongly with /i and vanish in the NDR scheme for 
fi = 5.5 GeV and = 6.3 GeV respectively. The corresponding values in the 
HV scheme are /i = 7.5 GeV and /i = 8.6 GeV . 

• There is no value of fi = fif in the full range considered for which £i(/i) and 
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Figure 2: (/i) in the NDR and HV schemes. 

£8 (a*) or equivalently £f F (/i) and £2 (//) simultaneously vanish. We also observe 
contrary to expectations in that is not necessarily smaller than e 8 ([i). 

In fact the large iV arguments presented in |7j that £i(p) = (9(l/iV 2 ) and 
£ 8 (/i) = 0(1 /N), imply strictly speaking only that the /i-dependence of £s(f-i) 
is much stronger than that of Ei(fi), which we indeed see in figs. |I] and |2|. The 
hierarchy of their actual values is a dynamical question. Even if the large Ac- 
counting rules = 0{1/N 2 ) and Eg(fi) = 0(1/ N) are true independently of 
the factorization hypothesis p0[ plf| , it follows from our analysis that once the 
generalized factorization hypothesis is made, the extracted values of £, violate 
for some range of \x the large-N rule E\ <C e$. 

We can next investigate for which renormalization scheme characterized by k± the 
factorization is exact at fif = = 4.8 GeV. We call this choice the "factorization 
scheme" (FS). Using the central values in (|50| ) and A^^L = 225 MeV we find by means 
of (fH) and (g§) 

k+ = 13.5 , k_ = 3.9 (FS). (51) 

These values deviate considerably from the NDR values k± = and the HV values 
k± = =f4. Yet one can verify that for these values J + = 6.13 and J_ = 1.17 and 
consequently in this scheme the NLO corrections at fi = remain perturbative. In 
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Table 1: CiJitA as functions of \i for different schemes and = 225 MeV. 









/i[GeV] 


NDR 


HV 


FS 


NDR 


HV 


FS 


2.5 


0.046 


0.035 


-0.033 


0.102 


0.144 


0.075 


5.0 


0.065 


0.059 


0.001 


0.022 


0.055 


-0.004 


7.5 


0.071 


0.067 


0.014 


-0.016 


0.013 


-0.041 


10.0 


0.074 


0.071 


0.021 


-0.039 


-0.013 


-0.064 



table we give the values of ^ NF (/i) for the NDR, HV and FS schemes. 

The numerical analysis presented here used as input the central values for a^ ff given 
in (|50|). As stressed in particular in f22"f , the strong model dependence of the form 



factors and large experimental errors preclude at present a precise determination of 
these parameters. Consequently when these uncertainties are taken into account, the 
differences between various schemes are washed out to some extent. Yet the general 
features of the results obtained for other numerical values of the pair (a® ff ,a| ff ) are 
very similar to the ones presented here. 

5 Generalized Factorization and N eS 



As pointed sometime ago in |13], Q and recently discussed in || g, [|, it is always 
possible to calculate the scale and scheme dependence of the hadronic matrix elements 
in perturbation theory by simply calculating the matrix elements of the relevant 
operators between the quark states. Combining these scheme and scale dependent 
contributions with the Wilson coefficients Cj(/i) one obtains the effective coefficients 
(7? which are free from these dependences. If one neglects in addition final state 
interactions and other possible non-factorizable contributions the decay amplitudes 
can be generally written as follows 

A = (H cS ) = ^VcMCfW™ + Cf (0 2 > tree ] , (52) 

where (Oj) tree denote tree level matrix elements. The proposal in || |], |9| is to use 
fl5"2"D and to apply the idea of the factorization to the tree level matrix elements. In 
this approach then the effective parameters a\ % are given by (||) with N cS treated as 
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a phenomenological parameter which models those non-factorizable contributions to 
the hadronic matrix elements, which have not been included into Cf s . In particular 
it has been suggested in |J that the values for N cS extracted from the data on 
two-body non-leptonic decays should teach us about the pattern of non-factorizable 
contributions. 

In particular when calculating the effective coefficients C? s , the authors of |§ U 
have included a subset of contributions to the perturbative matrix elements, which 
is sufficient to cancel the scale and scheme dependence of the Wilson coefficients. 
Unfortunately the results of such calculations are generally gauge dependent and 
suffer from the dependence on the infrared regulator and generally on the assumptions 
about the external momenta. 



Let us discuss this point in detail, following [23]. The Green function of the 



renormalized operator O, for a given choice of the ultraviolet regularization (NDR or 
HV for example), a choice of the external momenta p and of the gauge parameter A, 
is given by 

with 

r = r NDR ' HV + \r\ (54) 

The matrices f NDR ' HV depend on the choice of the external momenta and on the 
ultraviolet regularization, while f A is regularization- and gauge- independent, but de- 
pends on the external momenta. It is clearly possible to define a renormalization 
scheme in which, for given external momenta and gauge parameter, Tq(p) = 1, or in 
other words (0) Pt \ = (0) tree (this corresponds to the RI scheme discussed in p3|). 
However, the definition of the renormalized operators will now depend on the choice 
of the gauge and of the external momenta. If one were able, for example by means of 
lattice QCD, to compute the matrix element of the operator using the same renor- 
malization prescription, the dependences on the gauge and on the external momenta 
would cancel between the Wilson coefficient and the matrix element. If, on the con- 
trary, the matrix elements are estimated using factorization, no trace is kept of the 
renormalization prescription and the final result is gauge and infrared dependent. 

In m |] scale- and scheme-independent effective Wilson coefficients have been 
obtained by adding to Cj(/x) the contributions coming from vertex-type quark matrix 
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elements, denoted by fy and -yy. In particular 



(55) 



where the index j runs through all contributing operators, also penguin operators 
considered in |^, ^, 0]. 

It is evident from the above discussion that fy depends not only on the external 
momenta, but also on the gauge chosen. For example, in p, |9| the following result 
for fy is quoted: 
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(56) 



This result is valid in the Landau gauge (A = 0); in an arbitrary gauge, with the 
same choice of external momenta used to obtain (|56D one would get 



fy = fy(A 



0) + Xr v , 



(57) 



with fy(A = 0) given in ( pB[ ) and 
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The expressions for the full 10 x 10 f matrices in the NDR and HV schemes and in the 



Feynman and Landau gauges are given in p3fl , for a different choice of the external 
momenta. The results for the Landau gauge are given in [13 . 

Equation ( |5"7| ) shows that the definition of the effective coefficients advocated in 
[!• H H i s gauge-dependent. In addition, it also depends on the choice of the external 
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momenta. This implies that the effective number of colors extracted in ^, [J is 
also gauge-dependent, and therefore it cannot have any physical meaning. 

The gauge dependences and infrared dependences discussed here are not new. 
They appear in any calculation of matrix elements of operators between quark states 
necessary in the process of matching of the full theory onto an effective theory. A 
particular example can be found in |H| where the full gauge dependence of the quark 
matrix element of the operator (sd)v-A{sd)v-A has been calculated. However, in the 
process of matching such unphysical dependences in the effective theory are cancelled 
by the corresponding contributions in the full theory so that the Wilson coefficients 
are free of such dependences. Similarly in the case of inclusive decays of heavy quarks, 
where the spectator model can be used, they are cancelled by gluon bremsstrahlung. 
In exclusive hadron decays there is no meaningful way to include such effects in a 
perturbative framework and one is left with the gauge and infrared dependences in 
question. 

6 Summary 

In this paper we have critically analyzed the hypothesis of the generalized factoriza- 
tion. While the parametrization of the data in terms of a set of effective parameters 
discussed in 0-0 may appear to be useful, we do not think that this approach 
offers convincing means to analyze the physics of non-factorizable contributions to 
non-leptonic decays. In particular: 

• The renormalization scheme dependence of the non-factorizable contributions 
to hadronic matrix elements precludes the determination of the factorization 
scale [if. 

• Consequently for any chosen value of \ij = 0(m&) it is possible to find a renor- 
malization scheme for which the non-perturbative parameters e^s used in J7| to 
characterize the size of non-factorizable contributions vanish. The same applies 
to (aO introduced in the present paper. 

• We point out that the recent extractions of the effective number of colours N eS 
from two-body non-leptonic B-decays, presented in || §, U, while fi and renor- 
malization scheme independent suffer from gauge dependences and infrared 
regulator dependences. 
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A further problem in the generalized factorization approach is given by the pres- 
ence in many channels of operators that contribute only through non-factorizable 
terms. These contributions cannot be incorporated in the definitions of S\ and e§, 
and a more general parametrization is needed |25j. A typical example is given by 
charming-penguin contributions to B — > Kir decays |[26|| . 

We hope that our analysis demonstrates clearly the need for an approach to non- 
leptonic decays which goes beyond the generalized factorization discussed recently 
in the literature. Some possibilities are offered by dynamical approaches like QCD 
sum rules as recently reviewed in ||27|| . However, even a phenomenological approach 
which does not suffer from the weak points of factorization discussed here, would be 
a step forward. We hope to present some ideas in this direction in a forthcoming 
publication p8 . 
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